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ABSTRACT 

We present sub-milliarcsecond resolution total intensity and linear polarization VLBI images of 3C273, 
using concurrent 43 and 86 GHz data taken with the Very Long Baseline Array in May 2002. The structure 
seen in the innermost jet suggest that we have fortuitously caught the jet in the act of changing direction. The 
polarization images confirm that the core is unpolarized (fractional polarization < 1%) at 86 GHz, but also 
show well ordered magnetic fields (m ~ 15%) in the inner jet, at a projected distance of 2.3 pc from the core. 
In this strongly polarized region, the rotation measure changes across the jet by ~ 4.2 x 10 rad m~ 2 over an 
angular width of about 0.3 milliarcseconds. If the lack of polarization in the core is also attributed to a Faraday 
screen, then a rotation measure dispersion > 5.2 x 10 4 rad m~ 2 must be present in or in front of that region. 
These are among the highest rotation measures reported so far in the nucleus of any active galaxy or quasar, and 
must occur outside (but probably close to) the radio emitting region. The transverse rotation measure gradient 
is in the same sense as that observed by Asada et al and by Zavala and Taylor at greater core distances. The 
magnitude of the transverse gradient decreases rapidly with distance down the jet, and appears to be variable. 
Subject headings: galaxies: active — galaxies: jets — galaxies: magnetic fields — polarization — quasars: indi- 
vidual (3C 273) 
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1. INTRODUCTION 

Multi-frequency Very Long Baseline Polarimetry (VLBP) 
observations are crucial for determining the magnetic field 
structure and the distribution of Faraday rotating material in 
the innermost few parsecs of the cores of active galaxies and 
quasars. 

The fir st successful VL BP observations at 86 GHz were re- 
ported bv lAttridg el (1200 1 ). who presented polarized images of 
the well-known quasars 3C273 ansd 3C279. Here we report 
concurrent 43 GHz (A6.9mm) and 86 GHz (A3.5mm) VLBP 
of the 3C273 that allow us to investigate the distribution of 
Faraday rotating material near the core of this object, probing 
linear scales as small as a few light years 4 . 

3C 273 (J1229+0203) is a good source on which to attempt 
86 GHz VLBP, as it is one of the brightest quasars in the sky 
at that frequency. Also, its low redshift of z = 0.158 gives 
very high linear resolution (2.73 pc/mas) at the source. It is 
also a well studied superluminal source. For instance, Homan 
et al (2001) measured the proper motions of five jet compo- 
nents over six epochs in 1996 at both 15 and 22 GHz. Results 
ranged from 0.77 to 1.15 mas yr" 1 , corresponding to superlu- 
minal motions of (3 app = 7.9 - 1 1 .9. 

2. OBSERVATIONS AND CALIBRATION 

Concurrent 43 and 86 GHz observations of 3C 273 were 
performed on 9 May 2002 (epoch 2002.35) using the Very 
Long Baseline Array (VLBA) 5 over 10.5 hours. The VLBA 
was equipped with all ten 43 GHz receivers and seven 86 GHz 
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receivers (Fort Davis, Kitt Peak, Los Alamos, Mauna Kea, 
North Liberty, Owens Valley, Pie Town), effectively yielding 
the same resolution at each frequency. With a CLEAN beam 
of about 0.5 x 0.25 mas, the linear resolution at the source is 
about 6x3 light years. 

The data were recorded with a total bandwidth of 64 MHz, 
and were correlated at the National Radio Astronomy Obser- 
vatory (NRAO) in Socorro, NM. Approximately 1 .5 hours to- 
tal integration was spent on 3C273 at 43 GHz, and 3.0 hours 
at 86 GHz. These integration times lead to theoretical RMS 
noise values of ~0.5 mJy beam" 1 at 43 GHz, and ~2 mJy 
beam" 1 at 86 GHz. The observations of 3C273 were inter- 
leaved with those of the quasar 3C 279, which was a second 
target source and served to confirm our calibration. Results 
for 3C 279 will be presented in a separate paper. 

The NRAO AIPS and Caltech DIFMAP packages were 
used for calibration and hybrid imaging. Details about 
VLBP calibration and imaging may be found in ICottonl 
i 19931) and Roberts, Wardle, & Brown (1994). Methods 
to account for the complexities of 86 GHz calibration not 
present in lower frequency data, and to check the vali dity of 
the 86 GHz solutions are described in LAttridgd J2001I) and 
lAttridge. Greve. & Krichbauml ( fT9"9 9b). 

The calibration of the instrumental polarization (D-terms) 
was performed for 3C273 and 3C279 separately, and re- 
sults were compared. At 43 GHz, the typical D-term values 
were found to be ~ 0.03, and at 86 GHz they were — ■ 0.11. 
The r.m.s. difference between the D terms determined us- 
ing 3C273 and those determined using 3C279 was 0.02 at 
86 GHz. This is a measure of the errors in the determination 
of the D terms. Following the prescription in Appendix 1 of 
Roberts et al (1994), we estimate that the corresponding er- 
ror in fractional polarization in the images is ~ 0.003 (0.3%). 
This is consistent with the observed fractional polarization of 
the core component of 3C 273, which is < 1.0%. 

3. IMAGES 

Naturally weighted images of 3C 273 in total intensity, po- 
larized intensity, and fractional linear polarization at 43 GHz 
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FIG. 1. — Naturally weighted images of 3C273, epoch 2002.35, made with the VLBA at 43 GHz (left) and 86 GHz (right). In each panel we display (a) 
Total intensity distribution, (b) Linear polarization distribution, with tick marks showing the orientation \ °f me electric vectors, and (c) Fractional polarization 
distribution in grayscale (in color on astro-ph), over every other / contour. The restoring beams are shown in the lower left corners. At 43 GHz, the contours start 
at 15 mjy/beam and increase by factors of v2 up to peak values of 4547 mjy/beam (total intensity) and 630 mJy/beam (polarized intensity). At 86 GHz, the 
contours start at 25 mJy/beam and increase by factors of v2 up to peak values of 1815 mjy/beam (total intensity) and 256 mjy/beam (polarized intensity). 



and 86 GHz are shown in Figure 1 . Tick marks represent the 
orientation of the electric vectors of the polarized radiation. 
The orientations of the EVPAs at 43 GHz were derived using 
data from the "VLA/V LBA Polarization Calibration Page" 
JTavlor & Mversll200(3) . in which 3C273 was observed only 
one day prior to our own observations. 

Currently, there is no source to use as a calibrator to set the 
zero-point of the EVPAs at 86 GHz. We therefore arbitrarily 
chose to rotate tick marks in the 86 GHz P image so that the 
EVPA of component 7 is the same at 43 and 86 GHz. To 
see the EVPAs more easily, Figure 2 shows a zoom of the 
polarized emission from the core region at each frequency. 

Model fits to the u-v data are presented in Table 1; com- 
ponent labels for 86 GHz were matched to corresponding 
43 GHz components for ease of comparison. 

The 43 GHz total intensity image has a SNR of 2400, and 
a peak to lowest contour dynamic range of ~300:1. The 
86 GHz / image presented in Figure 2a has a SNR of 590, 
and a peak to lowest contour dynamic range of ~75:1. The 
corresponding numbers for the P images are about 7 times 
smaller. 

3.1. Total intensity images 

The 43 GHz images in Figure 1 show the well known core- 
jet structure of 3C 273. The core is component D and the jet 
stretches to the south-west (and continues for another 20,000 
mas in all wavebands). At this resolution, the jet is quickly 
resolved. At 86 GHz, the resolution is similar (because of the 
smaller array) but the lower sensitivity means that essentially 
only the first milliarcsecond of the jet can be seen. 



The component we label Ql/Wl is the same as component 
U4/K4 in Ojha et al (2004), who measured a proper motion 
of 9.9 ± 0.4c. Notice that the inner jet no longer points di- 
rectly towards this component. This is because the direction 
of the inner jet varies on fairly short timescales over a range 
of structural position angles (SPAs) between at least -110° 
and -140° (e.g. Krichbaum et al. 1990). Abraham & Romero 
(1999) proposed that the variation in orientation of the inner 
jet was periodic with a period of about 20 years, and they fit- 
ted a precessing jet model to the extant VLBI observations. It 
is now clear (Ojha et al 2003; Chen et al 2005) that the jet ori- 
entation varies on a timescale of months rather than decades, 
and is not periodic. 

The bright inner jet components in Figure 1, components 6 
and 7, are almost equidistant from the core at r ~ 0.9 mas, but 
at distinctly different structural position angles (SPAs). The 
northern component (6) is somewhat better aligned with the 
downstream jet components, and has a SPA of -123°. The 
southern component (7) is stronger, has a flatter spectrum, and 
a SPA of -144°. It is interesting that both are present simul- 
taneously. We are not aware of earlier images that show such 
a morphology, and we conjecture that we have caught 3C 273 
"in the act" of changing its jet direction. If this is so, compo- 
nent 7, because of its location and spectrum, is presumably the 
younger component, and is defining a new jet direction. Since 
it has already drawn level with component 6, we expect com- 
ponent 7 to be moving faster. Proper motion measurements 
will easily support or refute these conjectures. 

3.2. Polarization images 



43 and 86 GHz polarimetry of 3C273 
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TABLE 1 

Epoch 2002.35 Modelfit Results for 3C 273 
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NOTE. — Columns are as follows: (1) component name; (2) distance from eastern- 
most feature D; (3) position angle of separation from D; (4) total intensity; (5) fractional 
linear polarization; (6) orientation of the linear polarization position angle; (7) major 
axis of the model component; (8) minor axis of the model component; (9) orientation of 
the major axis. 

* All EVPAs at 86 GHz have been arbitrarily rotated so that the 86 GHz and 43 GHz 
values for this component are aligned. 



The polarized images shown in Figures 1 and 2 show that 
the inner jet is strongly polarized. Components Q6/W6 and 
Q7/W7 both exhibit fractional polarizations near 15% at both 
frequencies. Two years earlier, Attridge (2001) found a com- 
ponent at a similar location, whose fractional polarization was 
11% at 86 GHz. By contrast, linear polarization is not de- 
tected in the core (component D) at either frequency. The 
limits are 0.2% at 43 GHz, and 1% at 86 GHz. This also 
agrees with Attridge's (2001) observations. 

At lower frequencies, 3C273 also displays moderate lin- 
ear polarization in the jet, and little or no polarization in the 
core. These are properties common to many AGN (e.g. Taylor 
1998; Lister & Smith 2000; Pollack et al. 2003; Ojha et al. 
2004). 

4. DIFFERENTIAL FARADAY ROTATION ACROSS THE JET 

Although we are unable to determine the zero point for the 
EVPAs at 86 GHz, Figure 3 shows clearly that there is a large 
difference in Faraday rotation between jet components 6 and 
7. At 43 GHz the difference in EVPAs is xi~Xb = +32.7°. 
At 86 GHz the difference in EVPAs is xi~X6 = -33.9°. The 
difference between these numbers is 66.6° = 1 . 16 radians, and 
this is independent of the zero point at 86 GHz. There is there- 
fore an absolute difference in rotation measure between com- 
ponents 6 and 7 of ~ 3.2 x 10 4 rad m" 2 . In the source frame 
we multiply by (1 +z) 2 , which gives ~ 4.3 x 10 4 rad m~ 2 . 

The smallest possible rotation measures consistent with 
these data are -2.1 x 10 4 rad m" 2 for component 7 and 
+2.1 x 10 4 rad m" 2 for component 6. These are among the 
largest rotation measures ever observed. One or both of them 



are actually lower limits, because we have only determined 
their difference. The distance between components 6 and 7 is 
0.3 mas. If we express the rotation measure difference as an 
angular gradient, then it is 1 .3 x 10 5 rad m" 2 mas" 1 . 

Both Asada et al (2002) and recently Zavala & Taylor 
(2005) have also measured differences in rotation measure 
across the jet of 3C273 at the same location, about 7 mas 
from the core, but at different epochs (1995.9 and two epochs 
in 2000 respectively). Asada et al found a gradient of about 
80 rad m" 2 mas" 1 . Zavala & Taylor find a gradient of 500 rad 
m" 2 mas" 1 , and suggest that Asada et al measured a smaller 
value because of lower resolution (their observations were 
made at 5 and 8 GHz). The transverse gradient at 7 mas from 
the core has the same sign as we find at 0.9 mas (more positive 
rotation measure on the north side of the jet), but it is smaller 
by 2-3 orders of magnitude. (We would require simultaneous 
observations to make the previous statement more precise.) 

It is interesting to compare the rotation measures found in 
the middle of each slice where the signal to noise ratio is high- 
est and the effect of a gradient is negligible (so long as there 
is little depolarization). Asada et al measure about 350 rad 
m~ 2 at the center of their slice in 1995. This is consistent 
with the observations at 15 and 22 GHz throughout the fol- 
lowing year reported by Ojha et al (2004). Zavala and Taylor 
measure about 800 rad m" 2 at the same location some 4 years 
later. We suggest that all these measurements are in fact cor- 
rect, and that the rotation measure distribution varies in time. 



(a) 



(b) 
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FIG. 2. — Zoomed linear polarization distributions from Figure 1 reveal 
more clearly the orientation of the EVPAs as discussed in §4; (a) is at 43 GHz 
and (b) is at 86 GHz, both showing every second contour. 



4 



Attridge et al 



This would also be consistent with the results of deep mul- 
tiwavelength observations made by T. Chen (private commu- 
nication), and with observations of the core region (see next 
section). 

The Faraday rotation we observe must be external to com- 
ponents 6 and 7. Internal Faraday rotation will inevitably de- 
polarize the source (sometimes called "front-back" depolar- 
ization) and there is no sign of this (see Table 1). In general 
the fractional polarization drops by a factor of two when the 
observed rotation reaches ~ 45° (Burn, 1966). Here, the least 
depolarization would be observed if the rotation measures of 
components 6 and 7 were ±2.1 x 10 4 rad m" 2 . Even then, the 
predicted fractional polarizations at 43 GHz would be only 
~ 10%, whereas the observed values are close to 15%. This 
rules out internal Faraday rotation by material that fills the 
jet, but Faraday rotation in a sheath or boundary layer is still 
a possibility. Inoue et al. (2003) reached the same conclu- 
sion for the gradient observed by Asada et al., and the same 
argument can be applied to Taylor & Zavala's observations. 

5. FARADAY DEPOLARIZATION IN THE CORE 

No linear polarization is detected from the core component 
D in these observations. The limits are 0.2% at 43 GHz, and 
1% at 86 GHz. The high rotation measures observed less than 
1 milliarcsecond downstream from component D suggest that 
it is depolarized by differential Faraday rotation. This could 
be internal to the jet, but we favor an external screen because 
it is a natural extension of the screen that causes the Faraday 
rotation in front of components 6 and 7. If the intrinsic po- 
larization of component D is mo, and the variance in rotation 
measure along different lines of sight to component D is o\ M , 
the observe d polarization is given by m( A) = mo exp(-2cr| M A 4 ) 
JBurrI l966>. We do not know the value of mo, but a plausible 
range is from 0.7% to 7.0%, since this is the range of polar- 
izations observed by Nartallo et al (1998) at 270 GHz using 
the JCMT The corresponding minimum values for o-rm are 
1.6-8. Ox 10 4 radm" 2 . (The first value is derived from the 43 
GHz measurement.) We cannot be certain that the 270 GHz 
measurements by Nartallo et al (1998) actually refer to our 
component D. But if they do, we can derive an upper limit on 
<trm from the requirement that the core is significantly polar- 
ized at 270 GHz: a RM < 10 6 rad m" 2 . 

Note that this dispersion requires that the scale length for 
rotation measure fluctuations is smaller than component D it- 
self. Farther out, at 0.9 mas, the fluctuation scale length is 
larger than the size of components 6 and 7 (or they would 
show depolarization at 43 GHz) but smaller than the distance 
between them. Similarly, at 7 mas the fluctuation scale length 
is larger than the restoring beam but smaller than the width of 
the jet. 



Zavala and Taylor (2001) have also observed large and spa- 
tially resolved Faraday rotation in the core region of 3C273, 
which they found to be time-variable . At epoch 2000.07, 
they found a rotation measure of -2.6 x 10 3 rad m~ 2 in the 
north-east end of this region, near our component D, and 
+1.9 x 10 3 rad m" 2 in the south-east end, near the location 
of our components 6 and 7. Their measurements (which were 
at 15 and 43 GHz) do not match ours in any obvious way, 
probably attesting to the variability of this source. 

We note that if 3C 273 were at a more typical redshift of, 
say, z — 1, then components 6 and 7 would appear blended 
with the core, and could not be separated by current VLBI 
arrays. Then the apparent "core" would be polarized, and it 
would be very difficult to measure rotation measures as large 
as those reported here. Thus 3C273 may not be unusual in 
these properties; it is merely closer. 

6. CONCLUSIONS 

The results presented here demonstrate the value of polar- 
ization sensitive VLBI at 86 GHz. In particular, it allows us 
to probe the jet and its environment even closer to the cen- 
tral engine, and on the smallest available linear scales. It can 
also reveal very high rotation measures that cannot be seen at 
longer wavelengths because of source or bandwidth depolar- 
ization. 

At 86 GHz, 3C 273 reveals two strongly polarized (- 16%) 
components within a milliarcsecond of the core that have ro- 
tation measures that differ by ~ 4.3 x 10 4 rad m" 2 . This gra- 
dient is transverse to the jet direction and it has the same sign 
but is 2-3 three orders of magnitude larger than the gradients 
found by Asada et al (2000) and by Zavala & Taylor (2005) at 
greater distances from the core. Extrapolating such a gradient 
to the core component readily accounts for its persistent lack 
of polarization even at 86 GHz. 

The lack of depolarization in the jet requires that the Fara- 
day rotation occurs in a region external to the visble jet itself, 
but presumably close to it. The data suggest that this Fara- 
day rotating region is variable in time as well as spatially, and 
may be a boundary layer that represents the interaction of the 
jet with its environment. Rotation measure observations at 
the highest frequencies therefore offer a new way of explor- 
ing this environment. 
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